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O. 

The quasi-two-dimensional organic semiconductor a-(BEDT-TTF)2i3 has an 
anisotropic linear dispersion with the zero energy gap near the Fermi level. Due to the 
vanishing density of states at the Fermi level, the Coulomb interaction is unscreened in 
this material. We study theoretically the effect of the long-range Coulomb interaction 
and the low energy/long- wavelength behavior of a- (BEDT-TTF) 2I3 using the renormal- 
ization group analysis. The nearly logarithmic enhancement of the velocity reshapes the 
tilted Dirac cones, and changes the low temperature behavior. We also show the the- 
oretical calculation for the site-selective spin susceptibility, which can be measured by 
an NMR experiment. 
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Salts of organic molecule BEDT-TTF show various electronic phases, such as Mott 
insulator, charge transfer insulator, semimetal, and superconductor. In addition, there 
are many types of crystal structures. a-(BEDT-TTF)2i3 is one of such organic com- 
pounds, 2-4 which consists of alternately stacked BEDT-TTF molecule layers and tri- 
iodide (I J ) anion layers. BEDT-TTF molecular planes constitute conducting layers, on 
which quasi-two-dimensional electronic system is formed. The symmetry of the crystal 
structure of a-(BEDT-TTF)2i3 is low with only the spatial inversion. The unit cell con- 
tains four BEDT-TTF molecules, three of which are crystallographically inequivalent 
[Fig. 1 (a)]. 

The band calculation predicts that a-(BEDT-TTF)2i3 has semimetalic Fermi sur- 
face at room temperature and ambient pressure. 1 The first order phase transition to 
charge-ordered insulating phase 5 occurs at about 135 K with drastic changes in suscep- 
tibility 6 and resistivity. 7 At higher pressure, this charge-ordered insulating phase will 
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Fig. 1. (Color online) (a) The lattice structure of the conducting BEDT-TTF layer. A, B, and C 
denote the inequivalent sites. The inversion centers are also depicted, (b) The effective energy dispersion 
near the tilted Dirac cone (£ = +1). The Dirac cone is largely tilted in the x direction. The units in 
the graph are A -1 for momentum and eV for energy. 



gradually be suppressed, and completely vanishes at 1.5 GPa. It has been revealed 
that there exists a state with anisotropic linear dispersion near the Fermi level by the 
band calculation. 8-10 According to this calculation, a-(BEDT-TTF) 2 i3 has two strongly 
tilted Dirac cones, and the tilting is caused by the next-nearest-neighbor hopping. Ex- 
perimental results, e.g. T 2 dependence of the carrier density, 7 ' 11 are consistent with this 
linear dispersion. 

Graphene 12,13 is a representative material with Dirac cones. Unlike a- (BEDT- 
TTF) 2I3, graphene has a purely two dimensional electron system with isotropic Dirac 
cones. The band crossing point is just on the Fermi level, and vanishing density of states 
leads to the unscreened long-range Coulomb interaction. It has been analyzed by the 
renormalization group (RG) approach, and the Dirac cone reshaping due to the renor- 
malized Fermi velocity is observed by Shubnikov-de Haas oscillations as the change in 
the cyclotron frequency. 14 

Compared to graphene, the narrow band width of a-(BEDT-TTF)2i3 strengthen 
the electron correlation effect. It leads to the enhancement of the Dirac cone reshaping, 
and its effects are expected to be observed more easily compared to graphene. In this 
paper, we analyze the long-range Coulomb interaction effect on the system with tilted 
Dirac cones using the RG approach, and calculate the spin susceptibility as a physically 
measurable quantity. 

We start with the following generalized Weyl Hamiltonian describing the tilted mass- 
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less Dirac cones: 15 

U(k) = £W ■ k + V X k X <J X + VykyCTy, (1) 

where £ = ±1 denotes the valley degeneracy and we put h — 1. For the moment, we 
consider the £ = +1 valley. The parameter w determines the tilt of the anisotropic 
Dirac cone. The energy of this model [Fig. 1 (b)] is 

E ± (k) = wk±y/v*k* + v*kl (2) 

We assume that the parameters satisfy the relation 

WA \(<] 2 <1. (3) 



V x I \Vy 



This condition ensures that the system has a point node. 

We consider the anisotropic long-range Coulomb interaction 

V( Q ) = -gL^ (4) 

as a perturbation to the system. The anisotropy is reflected in the factor 77, and e is the 
dielectric constant. The unperturbed Green's function is 

G (k, u) = \ — . (5) 

U — W ■ K — V X K x a x — VyKyOy 

The RG analysis for two-dimensional systems often treats Coulomb interaction with 
large- iV expansion. 16,17 In the method for the isotropic system, the Coulomb propagator 
D (k, oj) is modified by adding one-loop fermion bubble diagram with TV fermion species: 

(2|*| + — 7 ==) . (6) 

The dressed term has importance in the strong coupling limit, but in the weak coupling 
case, it gives only small correction to the result. We concentrate our analysis on the 
low temperature region, where the running coupling constant becomes smaller than the 
bare value, so the dressed term is neglected in the following analysis. 

With the RG analysis, the parameters v and w are modified by the electron self- 
energy T,(k,uj). We calculate the self-energy to one-loop order. The one-loop order 
self-energy EW(fc,w) [Fig. 2] is evaluated by 

W(k,u) = i J ^-0- 2 G o (p,co + co')V(k-p). (7) 

The momentum integral is done in the momentum shell Ae~ l < \p\ < A. After some 
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Fig. 3. (Color online) The numerical solutions to the RG equations. Both v x and v y show almost 
completely logarithmic dependence to the momentum scale. 



calculation, we obtain 

= -£V X lk X <J X + ^VylkyCTy, (8) 

where the coupling constants a x and a y are defined by 

a - 2 ^ ^ C ° s2 dd ° (9a) 

£7T Jo (v I cos 2 9 + v 2 sin 2 9) 2 (cos 2 9 + rj sin 2 8) 2 

2e 2 p sin 2 8dd 

a y = — / — j — 3-. (9b) 

£n Jo (v 2 cos 2 9 + v 2 sin 9) 2 (cos 2 6* + 77 sin 0)5 

Then, the RG equations for v x and f y are 

£ =>, (10a) 

*k => (10b) 
To one-loop level of the self-energy, the tilting parameter w is not renormalized, and 
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Fig. 4. (Color online) The energy distribution in the fc-plane. The cutoff wavelength circle A = 
0.8 A 1 is shown by the dashed line. The solid lines denote the constant energy curves, plotted every 
0.02 eV. 

stays constant. The discussions above are unchanged for the £ = — 1 valley. 

The numerical solutions to Eq. (10) are shown in Fig. 3. The initial values at the 
cutoff momentum are v x = 0.0515, v y = 0.0439, w x = -0.0389, w y = 0.0048 (in eV A), 18 
and we set s — 10 and r\ — 1. v x and v y show nearly logarithmic dependence on the 
momentum scale like the system with the isotropic Dirac cone. 

The site-selective spin susceptibility \a ( a — A, B, and C) is described by 

Xa = J" deD a (e) , (11) 

where D a (e) is the site-dependent density of states (DOS), and /(e) is the Fermi dis- 
tribution. 18 The site-dependent DOS has the angular dependence in fc-plane. The DOS 
of site C mainly comes from the gentle slope of the tilted Dirac cone (9 ~ 0), and 
the DOS of site B comes from the steep slope (9 ~ ir). Site A has almost no angular 
dependence. Thus, the site-dependent DOS for each site shows different angular depen- 
dence. Site C, which corresponds to the gentle slope of the tilted Dirac cone, gives the 
largest contribution to the spin susceptibility, because it has the largest density of states 
among the three sites. On the contrary, site B brings the smallest spin susceptibility. 
Although the site-dependent DOS leads to the different magnitude of the site-selective 
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Fig. 5. (Color online) The theoretical calculation of the site-selective spin susceptibility. The solid 
lines and dashed lines describe the spin susceptibility for the interacting and the noninteracting case, 
respectively. 



spin susceptibility, the DOS for each site is proportional to the energy e. We can easily 
see that Xa oc T in the noninteracting system. 

For a two-dimensional system with linear dispersion, the density of states (DOS) 
D(e) is proportional to the energy e. However, the enhancement of the velocities by the 
RG analysis changes the energy dispersion, and sufficiently below the cutoff energy, the 
density of states as a function of the energy is suppressed compared to the noninteracting 
one. This scheme is valid for a-(BEDT-TTF)2l3 within the temperature range where 
the dispersion is well approximated by a linear one. 

The theoretical result for the site-selective spin susceptibility is obtained by numer- 
ical calculation, and it is shown in Fig. 5. We set the cutoff wavelength A = 0.8 A . 
Compared to the noninteracting result, which shows linear dependence of Xa in T, the 
RG analysis brings the reduction of the spin susceptibility Xa at low temperatures. 
The validity of the linear dispersion approximation depends also on the angle 9. The 
temperature ranges where the linear dispersion approximation holds are T < 70 K for 
the gentle slope and T < 100 K for the steep slope. 

Another important behavior is that the characteristic temperature of the site- 
selective spin susceptibility Xa is different for each site. Here, the characteristic tem- 
perature indicates the point where the spin susceptibility rapidly grows. Although we 
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set the circular cutoff momentum, the energy at the cutoff depends on the momentum 
direction due to the tilting of the Dirac cone. This fact leads to the strong suppression 
of the site-B spin susceptibility. 

In summary, we have studied the effect of the long-range Coulomb interaction in 
the system with tilted Dirac cones using the perturbative RG analysis. The velocity 
enhances logarithmically as observed in the isotropic case like graphene. We calculate 
the site-selective spin susceptibility for the quasi-two-dimensional organic conductor a- 
(BEDT-TTF) 2 i3. The RG analysis reduces the site-selective spin susceptibility at low 
temperatures, and the characteristic temperatures are different for each site. 
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